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membrane flow are considered. The flow induces coherent collective precession of chiral molecules. 
Our theoretical study shows that splay interactions in this system lead to spatial redistribution of 
chiral molecules and formation of spiral waves and target patterns observed in experiments. 
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Studies of pattern formation in nonequilibrium soft 
matter are essential for understanding the operation of 
biological cells and for potential applications Bio- 
logical membranes including active molecular pumps or 
channels can develop shape oscillations and show persis- 
tent wave propagation |2, |3|. If a membrane contains 
rotating molecular motors, interactions between them 
may lead to the development of regular arrays with the 
crystalline order . Phase separation in two-component 
lipid layers is responsible for budding and replication of 
vesicles 0,0- Closely related to biomembranes, Lang- 
muir monolayers are formed by organic lipid or am- 
phiphilic molecules disposed on a liquid-gas interface Q • 
Nonequilibrium patterns of traveling orientation waves in 
illuminated two-component Langmuir monolayers, where 
illumination leads to transitions between different confor- 
mational states of molecules, have been experimentally 
and theoretically investigated || EJ 0, [uj. Recently, 
Tabe and Yokoyama have demonstrated that Lang- 
muir liquid-crystal monolayers including chiral molecules 
("molecular rotors") are easily brought to and main- 
tained at nonequilibrium conditions by transmembrane 
flows 12]. If there is a gradient of small molecules across 
a Langmuir monolayer, i.e. their concentrations in the 
liquid and the gas are different, this produces a flow of 
such molecules through the monomolecular layer. Such 
transmembrane flow gives rise to coherent collective pre- 
cession of molecular rotors. Experiments using reflected- 
light polarizing microscopy have revealed that the preces- 
sion is not uniform and complex orientational wave pat- 
terns are observed. This behavior is apparently universal; 
it has been verified for a number of chiral chemicals and 
different experimental conditions |l2|. 

In this Letter, we construct a phenomenological theory 
of spatiotemporal pattern formation in chiral Langmuir 
monolayers with transmembrane flows. An important 
role in such systems is played by splay coupling between 
local concentration of chiral molecules and the orienta- 
tional field E El E3- We show that, in the pres- 
ence of transmembrane flow, this coupling gives rise to 



nonequilibrium wave patterns in the orientational field 
and spatial redistribution of chiral molecules inside the 
monolayer. The target patterns, seen in the experiments 
0, should thus be accompanied by aggregation of chi- 
ral molecules in the periphery of the patterns. Other 
wave patterns, such as traveling stripes and rotating spi- 
ral waves, are also possible. 

We study a model of an orientationally ordered two- 
component Langmuir monolayer representing a mixture 
of chiral and achiral molecules (in the experiments ^3 
the chiral molecules were making up only 10% of the 
monolayer). The local state of the monolayer is de- 
scribed by the variable c, giving the local fraction of chiral 
molecules, and by the director vector n that represents 
the projection of the molecular tilt onto the monolayer 
plane. The Landau free energy of the system is 

F = J [^(Vn) 2 + k B Tclnc + k B T(l - c) ln(l - c) 

-f |G(Vc) 2 + AcV • n]dxdy. (1) 

The first term corresponds to the elastic energy of orien- 
tational ordering [K is the Frank elastic constant). The 
next two terms determine the lattice-gas entropy con- 
tribution to the free energy (T is the temperature and 
k B is the Boltzmann constant), and the following term 
(with the coefficient G) takes into account weak ener- 
getic interactions between chiral molecules which favor 
their uniform spatial distribution. The last term in the 
expression for free energy describes splay interactions in 
the system. It provides coupling between the scalar con- 
centration field c and the vector orientational field n 5 
the parameter A specifies the strength of splay interac- 
tions fl(| . 

The kinetic equation for the local concentration c of 
chiral molecules is 
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where D is their diffusion constant. The kinetic equations 
for the director field n are 
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In addition to the relaxation terms (T is the relaxation 
rate constant for orientational ordering), we have phe- 
nomenologically included into these equations, following 
Ref . [l^ . a term that describes planar precession of the 
director vector. This precession is caused by the trans- 
membrane flow and its frequency Q is linearly propor- 
tional to the flow intensity (as seen in the experiments 
Because of the flow terms, the system cannot relax 
to the state of thermal equilibrium and oscillations and 
active wave propagation become possible. 

Rescaling time and spatial coordinates as t ^ 
£(/c#Tr) _1 and r — > r (K/kBT) 1 ^ 2 and using the angle 
variable <\> defined by n = (cos 0, sin 0), kinetic equations 
(0 and © can be written in the form 

c = v [V 2 c - gV (c (1 - c) V 3 c) + AV (c(l - c) V(V • n)] , 

(4) 
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6 = V 2 6 - uj + A cos 6- sin <j>— . (5) 
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The coefficients in these equations are v = ^(i^r) -1 , g — 
GK~\ A = A/(k B TK)^ 2 and uj = ft(kBTT)- 1 . Note 
that the total amount of chiral molecules is conserved 
and average spatial concentration Co of these molecules 
is a parameter of the system. According to equation @. 
splay coupling to the director field leads to physical forces 
acting on chiral molecules and to the viscous flow of these 
molecules in the monolayer plane. On the other hand, 
spatial gradients of concentration c lead, according to 
equation J5J), to local rotation of the director vector n. 

When transmembrane flow is absent (uj = 0), these 
kinetic equations describe relaxation to thermal equi- 
librium. The stationary equilibrium state is uniform 
if splay interactions are sufficiently weak. If the splay 
interaction strength A exceeds the critical value A cr = 
[co(l — co)] -1 ^ 2 , the uniform state becomes however un- 
stable with respect to growth of spatial modes with the 
wavenumbers < k < & max , where /^ ax oc (A — X cr )/g. 
This instability has previously been investigated and is 
known to lead to the formation of an equilibrium peri- 
odic stripe pattern In this equilibrium pattern, both 
the local concentration and the director orientation are 
periodically varying along a certain direction. 

To investigate nonequilibrium pattern formation in- 
duced by transmembrane flow, numerical simulations of 
the model (0 and (0) using the explicit Euler scheme 
with constant coordinate and time steps were performed. 
In all simulations, periodic boundary conditions were ap- 
plied. 

When the transmembrane flow is introduced (uj ^ 0), 
we see that the equilibrium stationary stripe pattern be- 
gins to move at a velocity that increases with the flux 
intensity uj. Figure 1 shows profiles of concentration and 
azimuthal angle across a traveling stripe pattern. Note 
that the shape of the stripes and their spatial period are 
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FIG. 1: Profiles of concentration c (solid) and azimuthal 
angle 0/27T (dashed) in a traveling stripe pattern for Co = 
0.2, v = 0.1, g = 1, A = 3, and u — 0.005. The arrow shows 
the direction of motion. 



not significantly different from the respective equilibrium 
pattern at uj = 0.. 

If a simulation is started with random initial condi- 
tions for the azimuthal angle field </>, regular stripes are 
not formed. Instead, the system undergoes relaxation 
to an equilibrium state with many spiral-shaped orien- 
tational defects. Application of the transmembrane flow 
to a system in this state leads, after a transient, to com- 
plex self-organized wave patterns. Several examples of 
such patterns for different parameter values are shown in 
Fig |21 and movies |l8j . 

The central region in the pattern shown in Fig. 2a is 
periodically emitting orientation waves. Repeated gener- 
ation and outward propagation of these waves is seen in 
the space-time diagram Fig J3] which displays evolution of 
the azimuthal angle distribution along the diagonal cross 
section indicated by the dashed line in Fig. 2a. The waves 
spread out in the large central region and run into the 
periphery part of the pattern, occupied by stripes with a 
shorter wavelength. The corners of the medium in Fig. 2a 
are occupied by spiral-shaped stripe structures (because 
periodic boundary conditions are used, they represent 
four parts of the same compact pattern). This stripe 
structure represents a wave sink, as can be seen from the 
space-time diagram in Fig|3| (the sink occupies the upper 
part of this diagram) and by examining the respective 
movie While a target pattern is seen in the cen- 

ter for the azimuthal angle distributions (right panel), a 
spiral wave occupying the central region is seen in the 
concentration distribution (left panel in Fig. 2a). 

Figure 4a displays a superposition of three subsequent 
snapshots of concentration profiles of chiral molecules 
along the diagonal line in Fig. 2a. Additionally, we 
show here the concentration profile c(x) smoothed over 
the spatial scale Ax = 56 (roughly the stripe period) 
and averaged over the time interval of 10000. Although 
spatial variations are substantial, we find that, on the 
average, the central region of the target pattern is de- 




FIG. 2: Distributions of concentration (left panels) and az- 
imuthal angle (right panels) in self-organized wave patterns 
obtained starting from random initial conditions for systems 
with the parameters: (a) Co — 0.1, A — 3.4, a; — 0.01, v — 
10, g = 10 and (b) c = 0.1, A = 3, oo = 0.015 v = 0.01, g = 1 
and (c) c 0.9, A 3.4, cj 0.015., v 0.01, # 10, The 
linear size of the medium is L — 800 (a) and L — 200 (b and c). 
The concentration is displayed in gray scale with the darker 
color corresponding to lower concentration values. 




FIG. 3: Space-time diagrams displaying evolution of the 
azimuthal angle along the diagonal cross section shown by 
dashed line in Fig. 2a. Time runs from left to right in the hor- 
izontal direction, the total shown time interval is T — 10000. 
The target pattern and the wave sink (above) are seen. 
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FIG. 4: (a) Snapshots of concentration profiles c(x,t) along 
the diagonal dashed line in Fig. 2a for three subsequent time 
moments (color online: red, green, blue); the black solid curve 
shows the average smooth profile c(x). (b) Effective local crit- 
ical splay interaction A cr (x) along the same line. 



pleted of chiral molecules which become concentrated in 
the curled stripes in its periphery. 

The spatial redistribution of chiral molecules allows 
to qualitatively explain the emergence of target-shaped 
wave patterns. When a stripe pattern, caused by splay 
interactions, is formed, it tends to hinder the angular 
rotation forced by the transmembrane flow. Therefore, 
oscillations develop only in the areas free from the stripes. 
As we have seen above, the critical splay interaction 
strength, needed for the formation of stripes, depends on 

1/2 

the concentration as A cr (c) = [c(l — c)] . If the con- 
centration is not constant and smoothly varies in space, 
the critical conditions will be determined by the local con- 
centration. As seen from Fig. 4a, the concentration varies 
rather rapidly inside the stripe regions. If we smooth the 
concentration profile and substitute c(x) instead of c in 
the expression for the critical splay interaction strength, 
we obtain the effective dependence X cr (x) displayed in 
Fig. 4b. For comparison, the horizontal dash line shows 
the value of A in the respective simulation. Thus, in- 
side the central area we have A < A cr (x), which explains 
why stripes are absent here. The boundary where the 
stripes first develop is roughly determined by the condi- 
tion A < X cr (x). In the region filled with stripes, A cr (x) 
lies below A, so that the uniform state is unstable. Of 
course, this argument is only approximate because the 
stripes are not at equilibrium. However, the effect of 
transmembrane flow on the stripe pattern is relatively 
weak, as we have already noticed. 

The simulation shown in Fig. 2a has been performed 
assuming strong diffusion of chiral molecules (y — 10). 
When diffusion is weaker {y = 0.01 in Fig. 2b), deple- 
tion of chiral molecules in the central region and their 
accumulation inside the stripe structure in the periph- 
ery become strongly pronounced. Now, the characteristic 
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wavelength of the stripes in the periphery of the target 
pattern is much shorter than the wavelength of the equi- 
librium stripe pattern for c = cq. This is because the 
wavelength of the stripe pattern decreases when the dif- 
ference | A— X cr (x) \ becomes larger. The stripes move only 
very slowly in the periphery under these conditions. The 
large central region can contain rotating spiral waves, as 
seen in Fig. 2b. The splay intensity strength A = 3 in this 
simulation was below the critical strength A cr = 3.33 • • • 
corresponding to Co = 0.1. This means that the uniform 
state remains stable with respect to small perturbations. 
However, strong initial perturbations can still lead in this 
case to the formation of steady wave patterns. The spiral- 
shaped stripe pattern in Fig. 2b is similar to the equilib- 
rium spiral patterns Note that a number of topo- 
logical orientational defects are present in the pattern 
shown in Fig. 2a, but all of them, except one, belong to 
the dense stripe region in the periphery. 

So far, only patterns for low average concentrations of 
chiral molecules have been discussed. In contrast to this, 
the wave pattern shown in Fig. 2c corresponds to a high 
concentration of chiral molecules (co = 0.9). Analyzing 
this pattern, several significant differences are seen. The 
concentration of chiral molecules is now increased inside 
the uniform central region and decreased in the region 
occupied by the stripes. The difference in the spatial 
distribution of chiral molecules for Co = 0.9 can be ex- 
plained if we notice that A cr (c) = [c(l — c)] 1 depends 
non-monotonously on concentration c and increases with 
concentration when c > 0.5. Therefore, inside the central 
region the formation of stripes is prevented because the 
condition A < A cr (x) again holds, now because the chi- 
ral molecules have aggregated in this region, pushing the 
achiral component into the stripe-filled periphery region. 

Thus, we have shown that splay interactions, based 
on coupling between the orientational and concentration 
fields, determine principal properties of nonequilibrium 
wave patterns in chiral Langmuir monolayers subject to 
the transmembrane flow. The theory explains target- 
shaped and spiral wave patterns observed in the experi- 
ments |l2|, [l9j . It relates the formation of such patterns 
to nonequilibrium spatial redistribution of molecular ro- 
tors, tending to aggregate in the areas occupied by slowly 
traveling, densely packed stripes. 

Biomembranes are closely related to Langmuir mono- 
layers and we expect that similar results should hold, 
under appropriate conditions, also for the membranes in- 
cluding a fraction of chiral molecules. The transmem- 
brane flow in such systems is created by a gradient of con- 
centration of small molecules or ions that leak through 
the membrane. The leakage may bring the membrane to 



nonequilibrium conditions, giving rise to traveling waves 
and complex self-organized wave patterns. Importantly, 
chiral molecules (and, possibly, some passive inclusions) 
can then be transported and spatially redistributed in a 
membrane as a result of wave propagation. 
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